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Abstract. A lattice computation of the B-meson decay constant and the mass of the b-quark to leading 
order in the heavy quark effective theory is presented. The involved renormalization problems are solved 
non-perturbatively, and the continuum limit is taken. In the quenched approximation the results reported 
here already offer an interesting numerical precision, which will be further improved in the near future. 



1 Introduction 

Not least by the influence of the phenomenologically very 
interesting programme of current experiments to inves- 
tigate CP-violation in the B-system 0, the study of B- 
meson physics has become a vivid area of research. To in- 
terpret the experimental observations within (or beyond) 
the standard model, matrix elements between low-energy 
hadron states must be known. But since these QCD ma- 
trix elements live in the strongly coupled sector of the the- 
ory, they naturally call for a genuinely non-perturbative, 
'ab initio' approach for their determination: the lattice 
formulation of QCD, which enables a numerical compu- 
tation of its low-energy properties through Monte Carlo 
evaluation of the underlying Euclidean path integral. 

Lattice QCD calculations with b-quarks can valuably 
contribute to precision CKM-physics by (over-) constrain- 
ing the unitarity triangle and help to obtain other phe- 
nomenologically relevant predictions. Examples for exper- 
imentally inaccessible key parameters that are important 
here are the B-meson decay constant and the mass of the 
b-quark, which are subject of the present report. In study- 
ing B-physics on the lattice, however, we face some par- 
ticular problems. A first one already originates in the b- 
quark itself, the mass of which is much larger than the 
inverse lattice spacings, 1/a, affordable in simulations on 
present-day computers even in the quenched approxima- 
tion (1/mb « 1/(4 GeV) < a « 0.07 fm): huge discretiza- 
tion errors would render a realistic treatment of B-systems 
with a propagating b-quark on the lattice impossible. 

This motivates to recourse to effective theories. Theo- 
retically most attractive is the heavy quark effective the- 
ory (HQET) whose Lagrangian in lattice formulation 
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is, to first order in the inverse heavy quark mass 1/m, 
formally identical to the continuum one. As a similar ex- 

* Invited talk at the International Europtiysics Conference 
on High, Energy Physics, EPS-HEP2003, July 17 - 23, 2003, in 
Aachen, Germany. 



pansion holds for the matrix elements in question, lat- 
tice HQET constitutes a systematic expansion in terms of 
l/rrib for B-mesons at rest |2] that also has a continuum 
limit order by order in the 1/m-expansion. 

Mainly because of two reasons, it has not received 
much attention in the past though: 

1. The rapid growth of statistical errors as the time sep- 
aration of correlation function is made large. This un- 
wanted feature is already encountered in the lowest- 
order effective theory (static approximation) and lim- 
its a reliable extraction of masses and matrix elements. 

2. The number of parameters in the effective theory does 
not only increase with the order of the expansion, but 
they have also to be determined non-perturbatively, 
since otherwise — as a consequence of the mixings 
among operators of different dimensions allowed in the 
cutoff theory (e.g. of j^Vh^^^V'h with ^/'jjDoV'h) — one 
is always left with a perturbative remainder that di- 
verges as a ^ 0. Hence, these power-law divergences 
cause the continuum limit not to exist unless the the- 
ory is renormalized non-perturbatively 3_. 

Here I summarize recent progress in both directions, 
which refiects in two concrete applications in the combined 
static and quenched approximation. These are a determi- 
nation of the Bg-meson decay constant, where a correction 
due to the finite mass of the b-quark is estimated by in- 
terpolating between the static result and Fd^ 0IE1, and a 
fully non-perturbative computation of the b-quark's mass 
based on the idea of a non-perturbative matching of HQET 
and QCD in finite volume as proposed in Refs. [HMZI- 



2 The Bg-meson decay constant 

On our way towards a precision computation of Fb^ in 
quenched QCD we employ a two-step strategy. First, 
the decay constant is calculated in lowest order of HQET, 
and then it is combined with available results for the pseu- 
doscalar decay constant i^ps(TOps) in QCD around the 
charm quark mass region by interpolation in l/mps- 
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The pseudoscalar decay constant at finite mass is re- 
lated to the renormalization group invariant (RGI) matrix 
element of the static axial current, 
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where 
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- Bg, through: 
O (1/M) . (2) 



Here, M denotes the RGI mass of the heavy quark and 
^MS QCD /l-parameter in the MS scheme. The renor- 
malization factor ^RGi, turning any bare matrix element 
of Ag*^* into the RGI one, has been non-perturbatively 
determined in jHI- Cps accounts for the fact that in order 
to extract predictions for QCD from results computed in 
the effective theory, its matrix elements are to be related 
to those in QCD at finite quark mass values. In this sense 
Cps translates to the 'matching scheme' |9I8| . which is de- 
fined by the condition that matrix elements in the (static) 
effective theory, renormalized in this scheme and at scale 
jjL = m, equal those in QCD up to l/m-corrections. In 
leading order it is given via the large-mass asymptotics 



^RGI - lim [ln(M/^j^)] 



7o/(2bo) 



^ps V^ps , (3) 



and thanks to the recent 3-loop computation of the anoma- 
lous dimension of the static axial current , the function 
C-ps{M / Ay^) = -Fpsy^mps/^RGi is now known perturba- 
tively up to and including 5*(m)-corrections. A numerical 
evaluation as explained in "W is shown in Figure ^ where 
one can also infer that the remaining perturbative uncer- 
tainty has become very small. 
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Fig. 1. Conversion factor to the matching scheme that trans- 
lates a RGI matrix element of Ag*^* at infinite mass to the one 
at finite mass. Its uncertainty is estimated to be smaller than 
2% (half of the difference between the 2- and 3-loop results). 

2.1 RGI matrix element in the static theory 

As mentioned before, heavy-light correlation functions on 
the lattice, from which B-physics matrix elements such 
as the B-meson decay constant in question are obtained 
at large Euclidean time, are affected by large statistical 
errors in the static approximation. Their noise-to-signal 
ratio grows exponentially with the time separation, and 
in particular for the Eichten-Hill action [5], 
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this ratio roughly behaves as exp{a;o(£'stat ^ ?ti„)} [TT] 
with -Estat the bare ground state energy of a B-meson in 
the static theory, diverging linearly in the continuum limit. 

To overcome this difficulty, we introduced in Ref. 0] 
a few alternative discretizations of the static theory that 
retain the 0(a) improvement properties of the action 
but lead to a substantial reduction of the statistical fluctu- 
ations. These new static quark actions rely on changes of 
the parallel transporters U{x,fi) in the covariant deriva- 
tive Dotph{x) — a~^[tph{x) — U^x — aO, 0)^phix — aO) ] of 
the form U{x, 0) W{x, 0), where now W{x, 0) is a func- 
tion of the gauge fields in the immediate neighbourhood 
of x,x + aO. Its best version employs 'HYP-smearing' that 
takes for W{x, 0) the so-called HYP-link, which is a func- 
tion of the gauge links located within a hypercube |12j. 
Comparing the noise-to-signal ratios in Figure |21 one can 
see that around xq w 1.5 fm more than an order of magni- 
tude can be gained in this case w.r.t. to the Eichten-Hill 
action and, in addition, the statistical errors grow only 
slowly as a;o is increased. Even more importantly, we ob- 
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Fig. 2. Noise-to-signal ratio of a B-meson correlation function 
in static approximation at a ~ 0.08 fm for various actions 
Bullets refer to the original Eichten-Hill action, while triangles 
correspond to the alternative discretization with HYP-links. 

served [4j (see also 13 ) quite the same, small lattice ar- 
tifacts with the new discretizations. 

In our computational setup to determine the bare ma- 
trix element (Bsl^o'''*!^) entering eq. we use the Schro- 
dinger functional (SF) formulation of QCD with non-per- 
turbatively 0(a) improved Wilson actions in the gauge 
and light (i.e. relativistic) quark sectors. For technical de- 
tails and the exact definitions of the correlation functions 
we refer to Refs. [30]. Here we only record that ^3] 
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^stat 
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,(a:o-T/2)B,tat(ano) 
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modulo volume factors, where /^*^* is a proper SF cor- 
relation function of the (0(a) improved) static axial cur- 
rent with the quantum numbers of a B-meson and /i is 
a corresponding boundary-to-boundary correlator, which 
serves to cancel the renormalization factors of the bound- 
ary quark fields. Moreover, we implement wave functions 
at the boundaries of the SF-cylinder to construct an inter- 
polating B-meson field that suppresses unwanted contam- 
inations from excited B-meson states to the correlators. 

So far, the bare matrix element lO has been calculated 
for three lattice spacings (a w 0.1 fm, 0.08 fm and 0.07 fm). 
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and the regularization dependent part of the factor Zrgi i 
which according to must be attached to get the RGI 
matrix element, has been computed for the new actions 
as it was done for the Eichten-HiU action in Ref. |H]- The 
continuum extrapolation quadratic in the lattice spacing 
of our results stemming from the static action with HYP- 
links is displayed in the right part of FigureOl To illustrate 
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Fig. 3. RGI matrix element of Aq*"' based on unimproved 
bare data from |15ll(i| (left) and on simulations of the 0(a) im- 
proved theory with the new discretization j4ll3j (right). A con- 
tinuum extrapolation of the latter yields rQ^^<pB.Gi = 1.74(13). 

the gain in precision and control of the systematic errors, 
we confront our 0(a) improved results with an analysis 
of older unimproved Wilson data for the bare matrix ele- 
ment, reproduced from ^ in the left part of the figure. 

2.2 Extrapolation in the heavy quark mass 

To finally arrive at a value for Fb_^ , we combine ^rgi , refer- 
ring to the static limit and thus being independent of the 
heavy quark mass, with numbers of _Fps in the continuum 
limit at finite values of the quark mass, which have been 
collected in the context of the (quenched) computation 
of Fd^ of In incorporating the mass dependence 

predicted by HQET, we are then led to extrapolate 

rf/^ Fps y'mps / Cps {M I Ay^) from the charm region to the 

static estimate Tq^^^rgi by a linear fit in 1/ (rgmps). This 
interpolation is shown in Figure^ where the zigzag error 
bands around the relativistic data indicate a small sys- 
tematic effect that is due to the mass dependence of the 
discretization errors in the decay constant near the charm 
quark mass as discussed in Refs. OE]- While an extrap- 
olation in 1/ (romps) from the charm region without the 
constraint through the static approximation would look 
similar, it is obvious that the interpolation is much safer, 
since extrapolating to the (quite distant) Bg-meson scale 
depends significantly on the functional form assumed. 

Using toBb = 5.4 GeV, rg = 0.5 fm and the numerical 
perturbative value of the matching factor Cps(Mb/ylj^) 
translating to finite b-quark mass, we find from the inter- 
polation to l/(romBg) in Figure^as our present result j2| 



roi^B, = 0.52(3) ^ Fb, = 206(10) MeV. 



(6) 



This number includes all errors except for the quenched 
approximation; the (unavoidable) scale ambiguity intro- 
duced by it can be estimated to be about 12%. 
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Fig. 4. Interpolation in the inverse heavy-strange meson mass, 
mps, between the RGI matrix element of the static axial cur- 
rent and relativistic data around the charm quark mass 

3 The b-quark's mass 

The second of the aforementioned problems that so far 
hampered the use of HQET on the lattice is the occur- 
rence of power-law divergences in the lattice spacing. It 
already shows up in the static approximation and thereby 
affects, for instance, the computation of the mass of the 
b-quark in leading order of HQET. In this case the ki- 
netic and the mass terms in the static action mix under 
renormalization and give rise to a local mass counterterm 
5m cx 1/a, the self-energy of the static quark, which im- 
plies a linearly divergent truncation error if one relies on 
an only perturbative subtraction of this divergence. There- 
fore, past lattice computations in the framework of HQET 
could not reach the continuum limit |19[I20) . 

A strategy for a solution to this longstanding problem 
was developed in Ref. [Jj, which now offers the possibility 
to perform clean, non-perturbative calculations in HQET. 
It basically consists of three parts that I want to briefly de- 
scribe in the following by sketching a (still ongoing) com- 
putation of the b-quark's mass as example [7ll21j : 

1. Renormalization of the effective theory amounts to re- 
late the parameters of the HQET Lagrangian to those 
of QCD, a step usually called matching. In order to 
realize the matching in a non-perturbative way, one im- 
poses matching conditions of the form ^^'^^"^(Lqj M) = 
^QCDj'^^ j^^-j jj-^ ^ physically small volume of linear 
extent Lq = 0(0.2 fm), where ^"QET ^qcd 
suitably chosen observables in HQET and QCD to 
be calculated with the aid of numerical simulations. 
The finiteness of the matching volume ensures that 
lattice resolutions satisfying am\, <^ 1 are possible 
and the b-quark can be treated as standard relativis- 
tic fermion, while at the same time the energy scale 
I/Lq = 0(1 GeV) is still significantly below mb and 
HQET applies quantitatively. In determining the pa- 
rameters of the effective theory from those of QCD 
via such a non-perturbative matching in finite volume, 
the predictive power of QCD is transfered to HQET. 
Of course, owing to the very construction of the effec- 
tive theory, it is clear that these matching conditions 
must also carry a dependence on the heavy quark mass, 
which is most conveniently identified with the (scheme 
and scale independent) RGI mass, M (see e.g. (2^). 
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In the concrete case of the b-quark mass computation, 
definite choices for the quantities have to be made 
to formulate a sensible matching condition between the 
quark mass in the full, relativistic theory (QCD) and 
HQET. Those are r(L,Mb), denoting the energy of 
a state with the quantum numbers of a B-meson but 
defined in a small volume of extent L, and FstiitiL) 
as its counterpart in (leading order of) the effective 
theory. As detailed in [71|53], both can be expressed 
as logarithmic derivatives of appropriate finite- volume, 
heavy- and static-light correlation functions, respec- 
tively, and numerically evaluated with high precision. 

2. Next we need to establish a connection to a physi- 
cal situation, where observables of the infinite- volume 
theory such as masses or matrix elements are acces- 
sible at the end. The accompanying gap between the 
small volume with its fine lattice resolution, where the 
matching of HQET and QCD is done, on the one side 
and larger lattice spacings (and also larger volumes) 
on the other is bridged by a recursive finite-size scaling 
procedure inspired by : the volume to compute the 
quantity .Tstat {L) in HQET is iteratively enlarged until 
one reaches a volume of linear extent L = 0(1 fm) so 
that, at the same resolutions a/ L (i.e. at the same bare 
parameters) met there, large volumes with L « 2 fm 
— to accommodate physical observables in the infinite- 
volume theory — eventually become affordable. Also 
note that, apart from terms of 0(l/(LoMb)"+^) if con- 
sidering HQET up to order n, any dependence on the 
unphysical small-volume physics is gone now. 

3. Finally, a physical, dimensionful input is still missing. 
In the case at hand this means to link the energy /Istat , 
which turns into the B-meson's static binding energy, 
-Estat, as the volume grows, to the mass of the B-meson 
as the physical observable in large volume whose nu- 
merical value is taken from experiment. 

To join the foregoing three steps, we have to recall that 
energies in the effective theory differ from the correspond- 
ing ones in QCD by a linearly divergent mass shift mbaro, 
which has its origin in the mixing of ip^^Doiph with the 
lower-dimensional operator V'hV'h under renormalization 
— the central problem we started from. As a consequence 
of its universality (i.e. its independence of the state), how- 
ever, mbarc obeys at any fixed lattice spacing 

rriB = Estut + "Ibaro , (7) 

r(L,Mb) = r,tat(i) + ^barc ■ (8) 

Imposing eq. ((HJ for L — Lq as the non-perturbative 
matching condition in small volume implicitly determines 
the parameter mbarc and may hence be exploited to re- 
place it in eq. Q). Then, after adding and subtracting 
a term r'stat(i2) (where L2 ~ 2^Lo w 1 fm with lattice 
spacings commonly used in large- volume simulations), the 
resulting equation may be cast into the basic formula 

TTlB = i?stat - il5tat(i2) + rstat(i2) " rstat(io) (9) 



I ^ in HQET 

r(io,Mb) 

a ^ in QCD 



a ^ in HQET 

O {A^/A'h) 



(with A a typical low-energy QCD scale), where the terms 
are just arranged such that (the unknown) mbaro can- 
cels in the energy differences AE = Estat — Atat(i2) and 
Atat(i2) — Atat(io), and the continuum limit exists sep- 
arately for each of the pieces entering eq. 

The entire heavy quark mass dependence is contained 
in r{Lo,M), defined in QCD with a relativistic b-quark. 
This mass dependence has been non-perturbatively map- 
ped out in Ref . |23| , where as a particular ingredient of the 
numerical calculation, which demands to keep fixed the 
(dimensionless) RGI heavy quark mass while approaching 
the continuum limit, the knowledge of several renormaliza- 
tion factors and improvement coefficients relating the bare 
to the RGI quark mass is required. Although they had al- 
ready been determined 22 25 , it was desirable to improve 
their precision and to estimate them directly in the bare 
coupling range relevant for our application. They were 
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Fig. 5. Continuum limit extrapolations of Lor{Lo, M) in 
a relativistic QCD and small volume (of linear extent Lo ~ 
0.2 fm) for representative values oi z = LqM |23I . 

thus redetermined in |23j and, as exemplified in Figure El 
performing controlled continuum extrapolations provides 
to = limg/ i^^Q Lor{Lo, M) as function of z = LqM. In 
view of the b-quark mass may now be extracted from 
the interception point of lo{z) with the combination 

Wstat = Lo»7iB-io{^stat(i2) - Tgtat (^o)} - ^O'^E' . (10) 



4.5 - 



4 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

"stat 






^(u,z) 

1 1 1 1 1 1 1 1 


1 


1 1 1 1 1 1 1 1 1 1 



5 5.5 ° " 6 6.5 z= 

Fig. 6. Solution of @ for the dimensionless RGI b-quark 
mass, LoMh- {u indicates that during the computation one 
works at a certain fixed value of the renormalized SF coupling.) 

The associated graph is given by Figure |S1 where for 
the time being we restricted the analysis to unimproved 
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Wilson fermion data for ai?stat from the literature 15 , 
resulting in LqAE = 0.46(5) (cf. the l.li.s. of Figured) . We 
presently obtain in static and quenched approximation |7] 

roMb 16.12(29) ^ m^(mp) = 4.12(8) GeV (11) 

up to corrections of OiA^/Mb) = 0{A/{LoMb)). From 
the — yet preliminary — r.h.s. of Figure [3 Qll^ one 
infers that, once the computation of aEstat with the static 
action discussed above (which also has linear 0(a) lattice 
artifacts removed) is finished, a continuum limit of LqAE 
with a (by a factor « 3) smaller error is in sight and will 
substantially enhance the accuracy of the result Hll|) . 
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Fig. 7. Subtracted, dimensionless Bs-meson energy evaluated 
from bare Wilson fermion numbers of [I5j (left) and from sim- 
ulations with the alternative discretization of [11 1211 (right). 

4 Conclusions and outlook 

This status report on actual work of our collaboration 
makes evident that, by virtue of (mainly two) recent ad- 
vances, non-perturbative calculations using the lattice reg- 
ularized heavy quark effective theory have reached a new 
quality. One important ingredient is the use of a modified 
static action which, for the first time, enables to compute 
B-meson lattice correlation functions with good statisti- 
cal precision in the static approximation for xq > 1 fm. 
As demonstrated both for the Bg-meson decay constant 
and for the b-quark mass, this represents a considerable 
improvement and has a great impact on the achievable 
precision in B-physics computations employing HQET.^ 
The determination of also applies the other promis- 
ing development, a general strategy how to solve renor- 
malization problems in HQET entirely non-perturbatively, 
taking the continuum limit throughout all steps involved. 

In the quenched approximation, where all the presented 
results refer to, the quoted uncertainties can (and will) be 
further reduced. Moreover, in mterpolating between data 
obtained in QCD and in the static limit, our result for Fb^ 
is almost independent of any effective theory. 

Finally it is worth to emphasize the interesting poten- 
tial of these methods for systematic and straightforward 
(albeit technically ambitious) extensions. Since it is one 
of the benefits of the theoretical concepts addressed here 

^ We note in passing that the results reported here agree well 
with those by a different new method that uses extrapolations 
in the heavy quark mass of finite- volume effects in QCD |2(il27| . 



that describing the b-quark by an effective theory circum- 
vents the need for prohibitively large lattices (because it 
completely eliminates the mass scale of the b-quark), they 
will very likely allow to also go beyond the static approx- 
imation by inclusion of 1/ mb-corrections as well as to in- 
corporate dynamical fermions without major obstacles. 
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